Break-up of shells under explosion and impact 
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A theoretical and experimental study of the fragmentation of closed thin shells made of a disor- 
dered brittle material is presented. Experiments were performed on brown and white hen egg-shells 
under two different loading conditions: fragmentation due to an impact with a hard wall and explo- 
sion by a combustion mixture giving rise to power law fragment size distributions. For the theoretical 
investigations a three-dimensional discrete element model of shells is constructed. Molecular dy- 
namics simulations of the two loading cases resulted in power law fragment mass distributions in 
satisfactory agreement with experiments. Based on large scale simulations we give evidence that 
power law distributions arise due to an underlying phase transition which proved to be abrupt and 
continuous for explosion and impact, respectively. Our results demonstrate that the fragmentation 
of closed shells defines a universality class different from that of two- and three-dimensional bulk 
systems. 
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Closed shells made of solid materials are often used 
in every day life, industrial applications and engineer- 
ing practice as containers, pressure vessels or combustion 
chambers. From a structural point of view aircraft vehi- 
cles, launch vehicles like rockets and building blocks of a 
space station are also shell-like systems, and even certain 
types of modern buildings can be considered as shells. 
The egg-shell as nature's oldest container proved to be a 
reliable construction for protecting life. In most of the 
applications shell-like constructions operate under an in- 
ternal pressure much higher than the surrounding one. 
Hence, careful design and optimization of structural and 
material properties is required to ensure the stability and 
reliability of the system. Closed shells usually fail due to 
an excess internal load which can arise either as a re- 
sult of slowly driving the system above its stability limit 
during its usage or service time, or by a pressure pulse 
caused by an explosive shock inside the shell. Due to 
the widespread applications, the failure of shell systems 
is a very important scientific and technological problem 
which has also an enormous social impact due to the hu- 
man costs arising, for instance, in accidental events. 

Fragmentation, i.e. the breaking of particulate mate- 
rials into smaller pieces is abundant in nature and under- 
lies several industrial processes, which attracted a contin- 
uous interest in scientific and engineeringresearch over 
the past decades 0, II II II II IE 0, H 11 E E E 
E. IfI Hfjl ITr| . Fragmentation phenomena can be ob- 
served on a broad range of length scales ranging from 
the collisional evolution of asteroids and meteor impacts 
on the astrophysical scale, through geological phenom- 
ena and industrial applications on the intermediate scale 
down to the break-up of large molecules and heavy nu- 



clei on the atomic scale. In laboratory experiments on 
the fragmentation of solids, the energy input is usu- 
allyachieved by shooting a projectile into a solid block 
ELM H 0. IE 13 making an explosion inside the sample 
Lj or by the collision of macroscopic bodies (free fall 
impact) El E EH EE!. Due to the violent nature 
of the process, observations on fragmenting systems are 
often restricted to the final state, making the fragment 
size (volume, mass, charge, ...) to be the main character- 
istic quantity. The most striking observation on fragmen- 
tation is that the distribution of fragment sizes shows a 
power law behavior, independently on the way of impart- 
ing energy, relevant microscopic interactions and length 
scales involved, with an exponent depending only on the 
djn^s^ality_of the system 0, H H || H Q, H SEE 
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efforts focused on the validity region and the reason of 
the observed universality in 1, 2, and 3 dimensions. De- 
tailed studies revealed that universality prevails for large 
enough input energies when the system falls apart into 
small enough pieces UMM E 111 B E EE E E 
E E E E, E HElSa , however, at lower energies 
a systematic dependence of the exponent on the input 
energy was evidenced E, E . Recent investigations on 
the low energy limit of fragmentation suggest that the 
power law distribution of fra gme nt sizes arises due to an 
underlying critical point E E E E E 

Besides the industrial and social impact of the fail- 
ure of shell like systems, they are also of high scientific 
importance for the understanding of fragmentation phe- 
nomena. Former studies on fragmentation have focused 
on the behavior of bulk systems in one, two and three di- 
mensions under impact and explosive loading, however, 
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hardly any studies have been devoted to fragmentation 
of shells j3g ■ The peculiarity of the break-up of closed 
shells originates from the fact that the local structure 
is inherently two-dimensional, however, the dynamics of 
the systems, the motion of material elements, deforma- 
tion and stress states are three-dimensional which allows 
for a rich variety of failure modes [38| . 

In this paper we present a detailed experimental and 
theoretical study of the fragmentation of closed solid 
shells arising due to an excess load inside the shell. Ex- 
periments were performed on brown and white hen egg- 
shells under two different loading conditions: fragmen- 
tation due to an impact with a hard wall and explosion 
by a combustion mixture have been considered resulting 
in power law fragment size distributions. For simplicity, 
our theoretical study is restricted to spherical shells such 
that a three dimensional discrete element model of spher- 
ical shell systems was worked out. In molecular dynamics 
simulations of the two loading cases, power law fragment 
mass distributions were obtained in satisfactory agree- 
ment with experiments. Based on large scale simulations 
we give evidence that power law distributions arise due to 
an underlying phase transition which proved to be abrupt 
for explosion and continuous for impact. Analyzing the 
energetics of the explosion process in the two loading 
cases and the evolution of the fragment mass distribu- 
tions we demonstrate that the fragmentation of closed 
shells defines a new universality class different from that 
of two- and three-dimensional bulk systems. 



I. EXPERIMENTS 

Hen eggs provide an excellent possibility for the study 
of fragmentation of thin brittle shells of disordered mate- 
rials with the additional advantages of being cheap and 
easy to handle, making the patience of scientists the only 
limiting factor for the subsequent improvement of the ex- 
perimental results. Our experiments were performed on 
ordinary brown and white egg-shells. In the preparations, 
first two holes of regular circular shape were drilled on the 
bottom and top of the egg through which the content of 
the egg was blown-out. The inside was carefully washed 
and rinsed out several times and finally the empty shells 
were dried in a microwave oven to get rid of all moisture 
of the egg-shell. The drying process proved to be essen- 
tial to ensure that the cuticula, which cannot be blown 
out, competely looses its toughness. 

In the impact experiments intact egg-shells are cata- 
pulted onto the ground at a high speed using a simple 
setup of rubber bands. The experimental setup provided 
a relatively high energy impact without the possibility 
of varying the imparted energy. The eggs are shot di- 
rectly into a plastic bag touching the ground so that no 
fragments are lost for further evaluation. 

In the explosion experiment initially the egg-shell is 
flooded with hydrogen and hung vertically inside a plastic 
bag. The combustion reaction is initiated by igniting the 




FIG. 1: Time evolution of the explosion of an egg-shell, con- 
secutive snapshots taken by a high speed camera. The time 
difference between the snapshots is 0.001 sec. 




FIG. 2: Time series of the impact of an egg-shell with the 
hard ground. The consecutive snapshots were taken by a high 
speed camera of 1 kHz. 



escaping hydrogen on the top of the egg. The hydrogen 
immediately reacts with the Oxygen which is also drawn 
up into the egg through the bottom hole, mixing with 
the remaining hydrogen. When enough air has entered 
to form a combustible mixture inside the egg, the flame 
back-fires through the top hole and starts the very quick 
exothermic reaction. The experiment is carried out inside 
a soft plastic bag so that secondary fragmentations due 
to fragment- wall collisions do not occur, furthermore, no 
pieces were lost after explosion. Since the pressure which 
builds up during combustion can slightly be changed by 
the hole size, i.e the smaller the hole, the higher the 
pressure at the explosion, we performed several series of 
experiments with hole diameters d between 1.2 and 2.5 
millimeter. The limit values have practical reasons: a 
drilling nail of large diameter typically breaks the eggs- 
shell, on the other hand it is extremely difficult to blow 
out an egg through a hole of diameter 1 mm or less. 

It is possible to follow the time evolution of the ex- 
plosion and impact processes by means of a high speed 
camera under well controlled conditions. Three consecu- 
tive snapshots of the explosion process are presented in 
Fig. n taken by a camera of 1000 Hz frequency. The ig- 
nition took place at the top of the egg in Fig. ^t). The 
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FIG. 3: Comparison of fragment mass distributions obtained 
by explosion experiments with two hole sizes and the impact 
experiment to the simulation results. The inset shows a typ- 
ical scanned set of fragments. 



tions F(m) for impact and explosion experiments aver- 
aged over 10-20 egg-shells for each curve. For the impact 
experiment, a power law behavior of the fragment mass 
distribution 

F(m) ~ mr T (1) 

can be observed over three orders of magnitude where 
the value of the exponent can be determined with high 
precision to r = 1.35±0.02. Explosion experiments result 
also in a power law distribution of the same value of r 
for small fragments with a relatively broad cut-off for 
the large ones. Smaller hole diameter d in Fig. [3J i.e. 
higher pressure, gives rise to a larger number of fragments 
with a smaller cut-off mass and a faster decay of the 
distribution F(m) at the large fragments. Comparing the 
number of fragments obtained, the ratio of the pressure 
values in the explosions at hole diameters d — 1.2 and 
2.0 mm, presented in Fig. [3] was estimated to be about 
1.6. Note that the relatively small value of the exponent t 
can indicate a cleavage mechanism of shell fragmentation 
and is significantly different from the experimental and 
theoretical results on fragmenting two-dimensional bulk 
systems where 1.5 < r < 2 has been found H H, 
EH El El 13 IH IMH |H 0113, and from the three- 
dimensional ones where r > 2 is obtained |2Ll5ll6L l39ll40| . 



instant of back-firing and the initiation of combustion is 
captured in Fig.^)), while in Fig. QJ;) already the flying 
pieces can be seen. Based on the snapshots the total du- 
ration of an explosion is estimated to be of the order of 
1 millisecond. 

In the impact experiment the egg hits the ground in 
the direction of its longer axis, as it is illustrated by the 
picture series of Fig. [21 After hitting the ground (Fig. 

the egg suffers gradual collapse as it moves forward 
(Fig. yp — h) making the impact process relatively longer 
compared to the explosion. 

The resulted egg-shell pieces are then carefully col- 
lected and placed on the tray of a scanner without over- 
lap. In the scanned image fragments are seen as black 
spots on a white background and were further analyzed 
by a cluster searching code. In the inset of Fig. [31 an ex- 
ample of scanned pieces of an impact experiment is shown 
where the broad variation of sizes can also be noticed 
with the naked eye. A dusty phase of shattered pieces 
|42| was also observed in the experiments with fragment 
sizes falling in the order of the pixel size of the scan- 
ner. The mass m of fragments was determined as the 
number of pixels in the scanned image. Since shattered 
fragments were also comparable to normal dust pieces in 
the air, they were excluded in the analysis by setting the 
lower cut-off of fragment masses to a few pixels. 

As the main quantitative result of the experiments 
we evaluated the mass distribution of fragments F(m) 
which is defined so that F(m)Am provides the proba- 
bility of finding a fragment with mass falling between m 
and m + Am. Fig.[3Jpresents the fragment mass distribu- 



II. SIMULATIONS 

Most of the theoretical studies on fragmentation relay 
on large scale computer simulations since capabilities of 
analytic approaches are rather limited in this field due to 
the complexity of the break-up process. Over the past 
years the Discrete Element Method (DEM) proved to 
be a very efficient numerical technique for fragmentation 
phenomena ^illllllllMllllsilSIMEfllill since 
it has the ability to handle large deformations arising 
dynamically, and naturally captures the propagation and 
interaction of a large number of simultaneously growing 
cracks, which is essential for fragmentation. 

In order to investigate the fragmentation of spherical 
shells we constructed a three-dimensional discrete ele- 
ment model such that the surface of the unit sphere is 
discretized into randomly shaped triangles (Dclaunay tri- 
angulation) by thro wing po ints randomly and indepen- 
dently on the surface [4311441 ]. Based on the triangulation, 
the dual Voronoi tessellation of the surface is also carried 
out as is illustrated in Fig.0] The nodes of the triangula- 
tion represent point-like material elements in the model 
whose mass is defined by the area of the Voronoi polygon 
assigned to it [3(| l43l l44j . The bonds between nodes are 
assumed to be springs having linear elastic behavior up to 
failure. Disorder is introduced in the model solely by the 
randomness of the tessellation so that the mass assigned 
to the nodes, the length and cross-section of the springs 
are determined by the tessellation (quenched structural 
disorder). After prescribing the initial conditions of a 
specific fragmentation process, the time evolution of the 
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FIG. 4: Example of the Delaunay triangulation of the spher- 
ical surface. The dual Voronoi lattice is also shown in the 
lower left quadrant. 



system is followed by solving the equation of motion of 
nodes by a Predictor-Corrector method of fourth order 

m-fi = Ff + Ff xt + Ff , i = l,...N, (2) 

where Ff is the sum of forces exerted by the springs con- 
nected to node i, and Ff xt denotes the external driving 
force, which depends on the loading condition. To fa- 
cilitate the relaxation of the system at the end of the 
fragmentation process, a small viscous damping force Ff 
was also introduced in Eq. J2Jl. 

In order to account for crack formation in the model 
springs are assumed to break when their deformation e 
exceeds a certain breaking threshold e c . A fixed thresh- 
old value e c — 0.03 is set for all the springs resulting in 
a random sequence of breakings due to the disordered 
spring properties. The breaking criterion is evaluated 
at each iteration step and those springs which fulfill the 
condition are removed from the simulation. As a result 
of successive spring breakings cracks nucleate, grow and 
merge on the spherical surface which can give rise to a 
complete break-up of the shell into smaller pieces. 

Fragments of the shell are defined in the model as sets 
of nodes (material elements) connected by the remaining 
intact springs. The process is stopped when the system 
has attained a relaxed state, i.e. when there is no spring 
breaking over a large number of iteration steps. The 
main advantage of DEM is that it makes it possible to 
monitor a large number of microscopic physical quanti- 
ties during the course of the simulation which are hardly 



Parameter 


Symbol 


Unit 


Value 


Initial radius 


R 


m 


1 


Initial volume 


Vo 


m 3 


4.19 


Initial Surface 


A 


m 2 


12.56 


Shell thickness 


th 


m 


5 • 1(T 5 


Total mass 


Mtot 


kg 


0.816 


Num. of triangles 


N t 




« 44000 


Num. of nodes 


N n 




« 21000 


Mass density 


P 


kg/m 3 


1300 


Time step 


At 


s 


3 ■ 10~ 7 


Damping coeff. 


Id 


kg/s 


0.1 


Spring Young modul. 


Y 


N/m 2 


10 9 



TABLE I: Parameter values used in the simulations. 



accessible experimentally, providing a deep insight into 
the fragmentation process. With the present computer 
capacities, DEM models can be designed to be realistic 
so that the simulation results can even complement the 
experimental information extending our understanding. 
The most important parameter values used in our simu- 
lations are summarized in Table [I] 

In computer simulations two different ways of load- 
ing have been considered which model the experimental 
conditions and represent limiting cases of energy input 
rates: (i) pressure pulse and (ii) impact load. A pressure 
pulse in a shell is carried out by imposing a fixed inter- 
nal pressure P a from which the forces F^ xt acting on the 
triangular surface elements are calculated as 

F?** = PoAjfti, (3) 

where Aj denotes the actual area of triangle j and the 
force points in the direction of the local normal tij, see 
also Fig.0] The force F? xt is equally shared by the three 
nodes of the triangle for which the equation of motion Eq. 
@ is solved. Since the surface area of the shell increases, 
the expansion under constant pressure implies a contin- 
uous increase of the driving force and of the imparted 
energy. 

The impact loading realizes the limiting case of instan- 
taneous energy input by giving a fixed initial radially ori- 
ented velocity v to the material elements and following 
the resulted time evolution of the system by solving the 
equation of motion Eq. J3J). The control parameter of the 
system which determines the final outcome of the process 
are the fixed pressure P Q and the initial kinetic energy E a 
for the pressure pulse and impact loading, respectively. 

III. THE BREAK-UP PROCESS 

In the simulations, in both loading cases the spheri- 
cal shell is initially completely stress free with no energy 
stored in deformation. When a constant pressure is im- 
posed the total energy E tot of the shell increases due to 
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Av = AV/V 

FIG. 5: The kinetic Efa n , elastic E e i, and total Etot energies 
as a function of the relative volume change Aw. Open symbols 
stand for expansion under constant pressure while the filled 
ones characterize the impact loading. The sudden drop of the 
total and elastic energy at Av c indicates the rapid break-up 
of the system. For Etot of the pressure loading the thick solid 
line follows Eq. JSJj while for E e i the function given by Eq. 
@ was fitted with C = 312000 as a parameter. 




the work done by the filling gas 



Etot{V) = I P dV = P AV, 
Jv 



(4) 



where V denotes the actual volume during the expan- 
sion and AV is the volume change with respect to the 
initial state V . The total energy can be written as 
the sum of the kinetic energy of material elements Ekm 
and of the elastic energy E e i stored in deformation, 
Etot — Ek in +E e i, where E e i is proportional to the change 
AA of surface area A of the expanding sphere with re- 
spect to the initial surface A Q . Introducing the relative 
volume change Av — ^r- as an independent variable, 
the total energy and the elastic energy can be cast in the 
form 



Etot = PoVoAv, 
E e i = c\(Av + l) 



,1/3 



1 



(5) 
(6) 



where the surface change AA was expressed in terms 
of Av. Furthermore, the parameter C of the system 
depends on the properties of the triangulation and the 
characteristic physical quantities of springs (Young mod- 
ulus, length, thickness). It is interesting to note that 
there exists a specific pressure value P* below which the 
expansion always stops at a maximum volume change 



FIG. 6: Cracks on the shell surface. Final states of impact 
experiments at energies E /E c ~ 0.8 o) and 2.8 b). Time 
evolution of the cracking process under a constant pressure of 
Po/Pc ~ 4.0 c, d) until the final relaxed state is reached e), 
with a magnified view of fragments /). Particle positions are 
projected back to their initial state on the surface. Fragments 
are identified as shell pieces surrounded by cracks. 



A.v max depending on P , however, for P a > P* the ex- 
pansion keeps always accelerating. For a given P Q < P* 
the value of Av max can be determined from the condition 
Etot = E e i so that 



PoV Av„ 



C (Au„ 



1 



,1/3 



1 



(7) 



and P* can be identified as the highest pressure for which 
Eq. (7J can be solved for Av max . Usually Av max can only 
be realized at low pressure values, because at higher pres- 
sures the system suffers complete break-up much below 
Av maxi due to the finite strength of the springs. Fig. 
[5] illustrates the evolution of the total E to t, kinetic E^in, 
and elastic E e i energies as a function of Aw for both pres- 
sure and impact loading. In the case of pressure loading 
it can be observed that the total energy E tot extracted 
from the simulations agrees well with the analytic predic- 
tion of Eq. J5J . The numerical value of the multiplication 
factor C of the elastic energy was obtained by fitting the 
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expression Eq. JBJ) to the curve of E e i(Av) in the figure. 
Due to the constant pressure, the total force F acting 
on the shell is proportional to the actual surface area 
F ~ A ■ P Q so that the system is driven by an increas- 
ing force during the expansion process. Since the driving 
force F increases with a diminishing rate when approach- 
ing the limit volume change Av ma xi it follows that the 
pressure loading case is analogous to the stress controlled 
quasistatic loading of bulk specimens. According to the 
simulations, under pressure loading there exists a critical 
pressure P c below which the expansion always stops at 
a finite volume and the shell only suffers partial failure 
( damage ) in the form of cracks but it keeps its integrity. 
When the pressure exceeds P c , however, the system sur- 
passes the critical volume change Av c when abruptly a 
large amount of springs break resulting in the break-up 
of the system (fragmentation). Note that P c « P* . 

The critical volume change Av c where fragmentation 
sets in during the expansion can be identified by the lo- 
cation of the sudden drop of the elastic energy in Fig. 
El caused by the large amount of spring breaking which 
occurs in a very narrow Av interval, resulting in a rapid 
formation of cracks on the surface. The value of Av c 
is mainly determined by the fixed breaking threshold e c 
and the disordered spring properties. Since the shell is 
under constant pressure the nucleated microcracks can 
grow and join giving rise to planar pieces surrounded by 
a free crack surface (fragment), as is illustrated in Fig. 
Efc, d, e). First large fragments are formed which then 
break-up into smaller pieces until the surviving springs 
can sustain the remaining stress, see Fig. Et, d, e). For 
simplicity, in the simulations the pressure is kept con- 
stant even if the system has lost its integrity, which has 
formally the consequence that pieces of the shell formed 
in the final state of fragmentation keep accelerating under 
the action of a constant force which explains the increas- 
ing kinetic energy Ekm in Fig. El following fragmentation. 
The volume of the system is numerically calculated as the 
sum of the volume of pyramidal objects defined by the 
surface elements and the center of the sphere, which pro- 
vides a meaningful result even after break-up in Fig. El 
in the vicinity of Av c . The critical pressure P c , required 
to exceed the critical volume change Av c to achieve frag- 
mentation, can be estimated as P c = E e i(Av c )/ \V Av c ). 

When loading is imposed by an instantaneous energy 
input E a , there is no further energy supply, the total en- 
ergy of the system is either constant or decreases due to 
the viscous dissipation and the breaking of springs (see 
Fig. EJ) ■ Since the elastic energy E e \ is solely determined 
by the deformation, the curve of E e i and the critical vol- 
ume change Av c where break-up arises in Fig. El coincide 
with the corresponding values of the pressure loading. 
Similarly to the pressure loading case, simulations re- 
vealed that a critical value of the imparted energy E c can 
be identified below which the shell maintains its integrity 
suffering only damage, while exceeding E c gives rise to a 
complete fragmentation of the shell. The resulted frag- 
ments on the shell surface obtained in the fragmented 



regime can be seen in Fig.Et, b). 



IV. FRAGMENT MASSES 

To give a quantitative characterization of the break- 
up of shells and to reveal the nature of the transition 
between the damaged and fragmented states large scale 
simulations have been performed varying the control pa- 
rameters, i.e. the fixed pressure P , and the imparted 
energy E a over a broad range. The most important char- 
acteristics of our fragmenting shell system, that can be 
compared to the experimental findings is the variation 
of fragment masses when changing the control parame- 
ters. In the simulations two cut-offs arise for the fragment 
masses, where the lower one is defined by the single un- 
breakable material elements of the model and the upper 
one is due to the finite size of the system. 

For both types of loading above the critical point the 
typical fragment size obtained at the instant of break- 
up decreases with increasing control parameter, which 
can be described analytically in terms of an energy bal- 
ance argument similarly to the one given in Ref. |2fij |. 
The loading energy of a shell region of linear extension 
L and mass m ~ L 2 , i.e. the energy stored in the 
motion of particles separating the piece from its sur- 
rounding, can be written as [m/M to t] Eki n (Av c )L 2 = 
[Ekin(Av c ) /M to t] L 4 , where Eki n {Av c ) denotes the total 
kinetic energy of the shell at the instant of break-up and 
Mtot is the total mass of the shell. The separation of the 
piece from its surrounding costs energy proportional to 
the fragment surface ~ L. The equilibrium fragment size 
can be obtained by minimizing the sum of the loading 
and surface energy densities pe 



E kin (Av c ) 2 1 



-1/3 
kin 



(8) 



It has 



with respect to L, which results in L ~ E 
been shown in the previous section that at the criti- 
cal point P c , E c the total kinetic energy of the system 
when break-up occurs takes zero value E k i n {Av c ) = 0. 
It follows that above the critical point E k i n has a lin- 
ear dependence on the distance from the critical point 
so that E kin (Av c ) ~ (P - P c ) for P Q > P c , and 
E kin (Av c ) ~ (E - E c ) for E a > E c hold. Substitut- 
ing these results into Eq. (jHJl, the typical fragment mass 
at the instant of break-up can be cast into the form 

m ~ (Po-P c y 2/3 for P >Pc, (9) 
m ~ {E -E c y 2 ^ for E > E c . (10) 

Eqs. H9I1U|I express that the typical fragment mass ob- 
tained at the time of break-up decreases according to 
a power law with increasing distance from the critical 
point. The exponent of the power law is universal in the 
sense that it does not depend on specific material prop- 
erties of the shell. Later on during the fragmentation 
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FIG. 7: Average mass of the largest (M m ax/Mtot) and second 
largest (M^fx / Mtot) fragment normalized by the total mass 
as a function of imposed pressure P . The inset presents a log- 
log plot of {Mmax/Mtot) as a function of the distance from 
the critical point P — P c . In the main figure the value of P c 
is indicated by an arrow. 




FIG. 8: {Mmax/Mtot) and [M^/Mtot) as a function of the 
imparted energy. The inset presents a log-log plot of the 
largest mass as a function of \E a — E c \, where the increasing 
and decreasing branches characterize the damaged and frag- 
mented states, respectively. The location of E c is indicated 
in the main figure. 



process the elastic energy stored in deformation may re- 
sult in succesive breakings of the large fragments. Hence, 
it can be expected that Eqs. (|9I10|I describe the scaling 
behaviour of the largest fragments, which did not un- 
dergo substantial size reduction until reaching the final 
relaxed state. 

Largest fragments. To characterize the degree of frag- 
mentation, i.e. the size reduction achieved in the sim- 
ulations, we calculated the average mass of the largest 
(M max /M to t) and of the second largest (M^ a d x /M tot ) 
fragment normalized by the total mass as a function of 
the pressure P , and input energy E Q in the case of pres- 
sure and impact loading, respectively [34l l35|. The re- 
sults are presented in Figs.[7| and[H] It can be seen that 
in both cases the maximum fragment mass is a mono- 
tonically decreasing function of the control parameters 
P a and E Ql however, the functional forms are different 
in the two cases. Low pressure values in Fig. [7| result 
in a breaking of springs, however, hardly any fragments 
are formed except for single elements broken out of the 
shell along cracks. Hence, the mass of the largest frag- 
ment is practically equal to the total mass M to t of the 
system, while the second largest fragment is orders of 
magnitude smaller (damage). Increasing however the 
pressure above the threshold value P c the largest frag- 
ment mass becomes much smaller than the total mass, 
furthermore, in this regime there is only a slight differ- 
ence between the largest and second largest fragments, 
indicating the complete disintegration of the shell into 
pieces (fragmentation) . The value of the critical pressure 



P c needed to achieve fragmentation and the functional 
form of the curve of (M max /M to t) above P c was deter- 
mined such that (M max /M to t) was plotted as a function 
of the difference \P — P c \ varying P c until a straight line 
is obtained on a double logarithmic plot. The result is 
presented in the inset of Fig.0where a power law depen- 
dence of (M max /M t ot) is evidenced as a function of the 
distance from the critical point 

(M max /M tot ) ~ \P -P c \- a , for P >P C . (11) 

The exponent a = 0.66±0.02 was obtained in good agree- 
ment with the analytic prediction of Eq. @. Detailed 
studies in the vicinity of P c revealed a finite jump of 
both (M max /M tot ) and {M^/ M tot ) at P c which im- 
plies that fragmentation occurs as an abrupt transition 
at the critical point, see Fig.[7| 

In Fig.[S]the corresponding results are presented for the 
case of impact loading as a function of the total energy E a 
imparted to the system initially. The mass of the largest 
fragment is again a monotonically decreasing function 
of the control parameter, however, it is continuous in the 
entire energy range considered. Careful analyzes revealed 
the existence of two regimes with a continuous transition 
at a critical value of the imparted energy E c . In the 
inset of Fig.\^(M max /M to t) is shown as a function of the 
distance from the critical point \E Q — E c \ where E c was 
determined using the same technique as for P c . Contrary 
to the pressure loading, (M max /M tot ) exhibits a power 
law behavior on both sides of the critical point but with 
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FIG. 9: The average fragment mass as a function of the im- 
posed pressure P . The inset presents a log-log plot of the 
average mass as a function of the distance from the critical 
point P — Pa for pressure values P > P c . The value of P c is 
the same as in Fig.|7| 



FIG. 10: The average fragment mass (M2/M1) as a function 
of the imparted energy E . The two regimes can be clearly 
distinguished. The location of the critical point proved to be 
exactly the same as in Fig. |5] The inset shows a log-log plot 
of the average mass as a function of \E — E c \. 



different exponents 



experiment j 



(M max I Mtot) 
{M max /M tot ) 



\E a - E e f, 
\E n -E r \- ( 



for E < E c , (12) 
for E Q > E c . (13) 



The numerical values of the exponents were obtained as 
a = 0.66 ± 0.02 and j3 = 0.5 ± 0.02, above and below 
the critical point respectively. Note that the value of 
a coincides with the corresponding exponent of the pres- 
sure loading and is in a good agreement with the analytic 
prediction of Eq. (jlOJI . Below the critical point the sec- 
ond largest fragment is again orders of magnitude smaller 
than the largest one, which implies that in this energy 
range the shell suffers only damage in the form of cracks, 
while above the critical point the break-up of the entire 
shell results in comparable values of the largest and sec- 
ond largest fragment masses. At the transition point E c 
between the damaged and fragmented states the mass of 
the second largest fragment has a maximum, while the 
curve of the largest one exhibits a curvature change, see 

Fig. m 

Average fragment mass. More insight can be obtained 
into the fragmentation process by studying the so-called 
single-event moments of fragment masses 



(14) 



where Ml denotes the fcth moment of fragment masses m 
in the jth realization of a fragmentation process, n 3 (m) is 
the number of fragments of mass m in event j. The ratio 
of the second M 3 2 and the first M( moments provides 
a measure for the average fragment mass in a specific 



M J 



Mi 



(15) 



Averaging over simulations with different realizations of 

disorder the average fragment mass M — (^M^/M^j was 

obtained as a function of the control parameter of the 
system. 

Due to the abrupt nature of the transition from the 
damaged to the fragmented states at the critical pressure, 
under pressure loading M cannot be evaluated below P c . 
However, when P a exceeds the critical pressure P c the 
average fragment mass monotonically decreases in Fig. [5] 
The inset of Fig. [5] shows M as a function of the distance 
from the critical point \P Q — P c \ where the same value of 
P c was used as in Fig. [7| A power law dependence of M 
is evidenced as a function of \P„ — P r \ 



M ~\P -P C 



(16) 



for P Q > P c and the value of the exponent was obtained 
to be 7 = 0.8 ± 0.02. For impact loading M can be 
evaluated on both sides of the critical point with a sharp 
peak in the vicinity of E c which is typical for continuous 
phase transitions in finite systems, see Fig. 1101 A power 
law dependence of M on the distance from the critical 
point 



M ~ \E n - E r 



(17) 



is again revealed for E a > E c , which is illustrated in 
the inset of Fig. The value of the exponent was de- 
termined by fitting 7 = 0.79 ± 0.02, which practically 
coincides with the 7 value of pressure loading. 
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FIG. 11: Mass distribution of fragments at various energies 
below and above the critical point. 
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FIG. 12: Mass distribution of fragments at various pressure 
values. 



Fragment mass distributions. The most important 
characteristic quantity of our system which can also be 
compared to the experimental results is the mass dis- 
tribution of fragments F(m). Under impact loading for 
E a < E c we found that F(m) has a pronounced peak at 
large fragments indicating the presence of large damaged 
pieces, see Fig. ^2 Approaching the critical point E c the 
peak gradually disappears and the distribution asymp- 
totically becomes a power law at E c . We can observe in 
Fig. ^] that above the critical point the power law re- 
mains for small fragments followed by a cut-off for the 
large ones, which decreases with increasing E a . 

For pressure loading F(m) can only be evaluated above 
P c . The evolution of F(m) with increasing pressure is 
presented in Fig. 1121 where the mass distribution always 
shows a power law behavior for small fragments with a 
relatively broad cut-off for the large ones. For the pur- 
pose of comparison, a mass distribution F(m) obtained 
at an impact energy close to the critical point E c , and 
distributions at two different pressure values P of the 
ratio 1.6 are plotted in Fig. [3] along with the experimen- 
tal results. For impact an excellent agreement with the 
experimental and theoretical results is evidenced. For 
pressure loading, the functional form of F(m) has a nice 
qualitative agreement with the experimental findings on 
the explosion of eggs, furthermore, distributions at the 
same ratio of pressure values obtained by simulations 
and experiments show the same tendency of evolution, 
see Fig. 

Figs.ll3landll4ldemonstrate that by rescaling the mass 
distributions above the critical point by plotting F(m) ■ 

M as a function of m/M an excellent data collapse is 
obtained with 8 = 1.6 ± 0.03. The data collapse implies 



the validity of the scaling form 

F(m) ~ m T ■ /(m/M), (18) 

typical for critical phenomena. The cut-off function / 
has a simple exponential form exp (— m/M) for impact 
loading (see Fig. I13|) . and a more complex one containing 
also an exponential component for the pressure case (see 
Fig. I14|) . The average fragment mass M occurring in the 
scaling form Eq. (|18|l diverges according to a power law 
given by Eqs. I|16I17|) when approaching the critical point. 
The good quality of collapse and the functional form Eq. 
(|18fl also imply that the exponent r of the mass distribu- 
tion does not depend on the value of the pressure P a or 
the kinetic energy E Q contrary to the bulk fragmentation 
where an energy dependence of r was reported |28| . 

The rescaled plots make possible an accurate determi- 
nation of the exponent r, where r = 1.35 ± 0.03 and 
r = 1.55 ± 0.03 were obtained for impact and pressure 
loading, respectively. Hence, a good quantitative agree- 
ment of the theoretical and experimental values of the 
exponent r is evidenced for the impact loading of shells, 
however, for the case of pressure loading the numerically 
obtained exponent turned out to be somewhat higher 
than in the case of exploded eggs. 

V. DISCUSSION AND OUTLOOK 

We presented a detailed experimental and theoretical 
study of the break-up of closed shells arising due to a 
shock inside the shell. For the purpose of experiments 
brown and white hen egg-shells were carefully prepared 
to ensure a high degree of brittleness of the disordered 
shell material. The break-up of the shell was studied 
under two different loading conditions, i.e. explosion 
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FIG. 13: Rescaled plot of the mass distributions for imparted 
energies above the critical point E > E c . The dashed line 
shows the fitted power law with an exponent r = 1.35 ± 0.03. 




FIG. 14: Rescaled plot of the mass distributions of various 
pressure values above the critical point P > P c . The dashed 
line indicates the fitted power law with an exponent r = 1.55± 
0.03. 



caused by a combustible mixture and impact with the 
hard ground. As the main outcome of the experiments, 
the mass distribution of fragments proved to be a power 
law in both loading cases for small fragment sizes, how- 
ever, qualitative differences were obtained in the limit of 
large fragments for the shape of the cut-off. 

We worked out a discrete element model for the break- 
up of shells which provides an insight into the dynam- 
ics of the process by simultaneously monitoring several 
microscopic quantities in the framework of molecular dy- 
namics simulations. In the simulations two ways of load- 
ing have been considered, which mimic the experimental 
conditions and represent limiting cases of energy input 
rates: during an expansion under constant pressure P 
the shell is driven by an increasing force with a contin- 
uous increase of the imparted energy, while the impact 
loading realizes the instantaneous input of an energy E a . 
Simulations revealed that depending on the value of P Q 
and E a , the final outcome of the break-up process can be 
classified into two states, i.e. damaged and fragmented 
with a sharp transition in between at a critical value of 
the control parameters P c and E c . In the fragmented 
regime power law fragment mass distributions were ob- 
tained in satisfactory agreement with the experimental 
findings. Analyzing the behavior of the system in the 
vicinity of the critical point P c , E c , we showed that power 
law distributions arise in the break-up of shells due to an 
underlying phase transition between the damaged and 
fragmented states, which proved to be abrupt for explo- 
sion, and continuous for impact. 

Due to its unique characteristics, the break-up of shells 
defines a new universality class of fragmentation phenom- 
ena, different from that of the two- and three-dimensional 
bulk systems. Based on universality, our results should 



be applicable to describe the break-up of other closed 
shell systems composed of disordered brittle materials. 
Explosion of shell-like fuel containers, tanks, high pres- 
sure vessels often occur as accidental events in industry, 
or in space missions where also the explosion of complete 
satellites may occur creating a high amount of space de- 
bris orbiting about Earth. For the safety design of shell 
constructions, and for the tracking of space debris it is 
crucial to have a comprehensive understanding of the 
break-up of shells. Due to the universality of fragmen- 
tation phenomena, our results can be exploited for these 
purposes. 

In the fragmentation of bulk systems under appropri- 
ate conditions a so-called detachment effect is observed 
when a surface layer breaks off from the bulk and un- 
dergoes a separated fragmentation process [2^, H(| . This 
effect also shows up in the fragment mass distributions 
in the form of a power law regime of small fragments 
of an exponent smaller than for the large ones. Our re- 
sults on shell fragmentation can also provide a possible 
explanation of this kind of composite power laws of bulk 
fragmentation |28l l3r| . 
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